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ABSTRACT 

This paper details a numerical investigation of performance 
losses in low-thrust solar thermal rocket nozzles. The effects 
of nozzle geometry on three types of losses were studied; 
finite rate dissociation-recombination kinetic losses, two- 
dimensional axisymmetric divergence losses, and 
compressible viscous boundary layer losses. Short nozzle 
lengths and supersonic flow produce short residence times in 
the nozzle and a nearly frozen flow, resulting in large kinetic 
losses. Variations in geometry have a minimal effect on 
kinetic losses. Divergence losses are relatively small, and 
careful shaping of the nozzle can nearly eliminate them. The 
boundary layer in these small nozzles can grow to a major 
fraction of nozzle radius, and cause large losses. These losses 
are attributed to viscous drag on the nozzle walls and flow 
blockage by the boundary layer, especially in the throat 
region. Careful shaping of the nozzle can produce a 
significant reduction in viscous losses. 


1 INTRODUCTION 

Continued exploration and exploitation of space will require 
the development of more efficient and economical systems for 
access to and transportation through space. Due to its 
potential for high specific impulse (I sp ), solar thermal 
propulsion has emerged as a candidate to meet this role. A 
schematic of a generic solar thermal rocket is shown in Fig. 1. 

The system consists of mirrored concentrators which 
concentrate solar radiation on a metallic absorber. The 
absorber transfers the energy to a working fluid. This working 
fluid is heated and then thermally expanded through a 
contoured nozzle to produce thrust. Hydrogen, with its low 
molecular weight, produces high exhaust velocities and is 
commonly considered for the working fluid propellant. 
Practical system sizes and weights generally limit solar 
thermal design concepts to low thrust systems. 

These characteristics result in the hydrogen fueled solar 
thermal rocket having about twice the I sp of a traditional 
chemical system, but requiring much longer times to transfer 
payloads. Therefore, the primary mission envisioned for solar 
thermal propulsion is in an orbital transfer stage. This 
mission can be divided into two categories: 1.) orbital transfer 
between LEO and GEO; and 2.) planetary/asteroid flyby 
missions such as the Pluto fast flyby being developed by 



FIGURE 1 : SCHEMATIC OF A GENERIC SOLAR THERMAL 
TRANSFER VEHICLE. 


NASA. These two mission categories are differentiated by the 
requirements on AV, mass, and thrust of the system. The 
LEO/GEO transfer vehicle requires relatively low AV's (~ 
10,000 ft/s, or ~ 3,000 m/s), but high mass, and thrust on the 
order of 6 - 40 lbf (25 - 180 N) for Delta class vehicles, or 
approximately 2 lbf (9 N) for smaller vehicles such as 
Conestoga. The planetary flyby mission requires 

significantly larger AV’s, relatively low mass, and thrust 
levels on the order of 1 - 2 lbf (4 - 8 N). These categories are 
not definitive and can significantly overlap depending on 
specific mission requirements. 

The solar thermal system is conceptually simple, but there 
are significant design issues which must be addressed. These 
issues include the design of the concentrator, development of 
the absorber/heat exchanger, and the optimization of the 
nozzle to efficiently produce the required thrust. Research 
initiatives are gradually developing to examine each of these 
component issues. Based on preliminary system studies by 
the Air Force and industry, it appears that a significant drop in 
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I S p (10 - 15%) can occur in low thrust (1-10 lbf, ~ 4 - 40 N) 
systems (Shoji, 1993). The physical causes of this 
performance drop must be adequately understood in order to 
design a practical solar thermal propulsion system. 

A liquid hydrogen (LH 2 ) propellant can produce theoretical 
I S p values of 700-1100 seconds based on typical temperatures 
in the absorber/heat exchanger of 3500-6500°R (1944-3611 
K). The temperature is primarily limited by the melt 
temperature of the absorber material and the cooling 
techniques available. At these temperatures, the liquid 
hydrogen is vaporized and a fraction of the gas molecules 
dissociate. As this hot gas thermally expands through the 
nozzle, it cools and recombination occurs with an associated 
heat release. The process of dissociation, recombination and 
heat release can result in a net reduction of I sp (Sutton, 1992; 
Kim and Stubbs, 1993). The analysis of this behavior is 
complicated by the fact that this may also be a 
thermochemicaily nonequilibrium process. 

Other possible sources of I S p reduction can be attributed to 
viscous effects and flow divergence. The specific nozzle 
geometry can have a significant impact on these and the 
kinetic loss. Therefore, designing nozzles for optimum 
performance is based on careful tailoring of the geometric 
contours. This paper describes a fundamental study of the 
potential causes of I S p reduction in low thrust, solar thermal 
rocket nozzles. The main emphasis is an analytical/numerical 
study of how nozzle geometry impacts the viscous and 
thermochemical losses. The results are applicable to 
optimizing nozzle geometry for maximum performance. 


2 MODEL DESCRIPTION 


2.1 Numerical Model 

The primary tool used in this investigation was the Two- 
Dimensional Kinetics (TDK) code of Nickerson, et al. (1993). 
The code is a robust engineering level design and analysis tool 
which is used extensively in the propulsion industry. The 
TDK code provides various nozzle analysis options. The 
nozzle flow field can be solved as quasi-one-dimensional, or 
with an axisymmetric method of characteristics formulation. 
Frozen, nonequilibrium finite rate (kinetic), and equilibrium 
chemistry options are available. These solutions are referred 
to in the code and in this paper as follows: One-Dimensional 
Frozen (ODF), One-Dimensional Kinetic (ODK), One- 
Dimensional Equilibrium (ODE), Two-Dimensional Frozen 
(TDF), Two-Dimensional Kinetic (TDK), and Two- 
Dimensional Equilibrium (TDE). 

The TDK Boundary Layer Module (BLM) uses a first-order 
integral method to compute nozzle performance losses due to a 
compressible laminar and/or turbulent boundary layer. This 
approach consists of using the initial method of 
characteristics solution to calculate the distribution of the 
nozzle boundary layer displacement and momentum 
thicknesses. The nozzle surface is then displaced by the 
displacement thickness to produce a new effective nozzle 
shape which can be re-analyzed by the method of 
characteristics. The BLM technique includes axisymmetric and 
longitudinal curvature effects. In addition, the nozzle wall can 



be modeled as adiabatic or with a fixed distribution of wall 
temperatures or heat fluxes. 

TDK provides various options for defining 
convergent/divergent nozzle geometries including simple 
conical nozzles and bell nozzles. The bell nozzle design used 
in this study is based on a skewed parabola as shown in Fig. 2. 
The nozzle transitions from a constant area combustion 
chamber into a conical convergent section defined by a 
convergence angle 0 U . The throat is defined by two circular 
arcs, one upstream with radius R u , and one downstream with 
radius R^. The length of the upstream arc is controlled by R u 
and the convergence angle, 0 U . The length of the downstream 
arc is controlled by R^ and the divergence angle, 0^. The 
nozzle expansion is fixed by fitting a skewed parabola to the 
location and slope of the downstream end of the throat and the 
nozzle exit. All nozzle dimensions are normalized by the 
throat radius, R*. 


2.2 Nozzle Performance Analysis 

The specific impulse for a rocket is defined as 

F 

Isp = - (I) 

m p g 

where F is the thrust, m p is the propellant mass flow rate, and 

g is the gravitational constant. The specific impulse is 
therefore a measure of how efficiently a rocket produces thrust. 
As discussed by Sutton (1992) and O'Leary and Beck (1992), 
I S p loss mechanisms for typical bell nozzles can be divided 
into three primary categories: (1) chemical kinetics losses, (2) 
divergence losses, and (3) viscous losses. Using TDK 
analysis results, the effects of each of these mechanisms can 
be expressed as a specific impulse loss, AI g p. 

In a solar thermal nozzle, the liquid hydrogen propellant is 
heated in the absorber section. This thermal energy is used to 
vaporize the liquid, but also produces some fraction of gas 
molecule dissociation. It is desirable to convert a significant 
fraction of this chemical energy to gas kinetic energy at the 
nozzle exit by recombination. The ideal case is for the gas to 
attain chemical equilibrium at every point in the nozzle flow 
field. However, the expanding flow can accelerate gas 
molecules so rapidly that their residence time in the nozzle 
does not permit full equilibrium to be reached. This state of 
thermochemical nonequilibrium results in a loss in efficiency. 
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A measure of this effect can be calculated as the difference in 
the TDK one-dimensional equilibrium and one-dimensional 
kinetics (nonequilibrium) solutions. 

^sp(KIN) = I sp(ODE) " I sp(ODK) ^ 

Theoretically, the kinetics loss can be minimized by 
increasing the residence time of the fluid in regions where the 
temperature is decreasing and the pressure is still high enough 
to promote recombination. 

Divergence losses occur when portions of the exit 
flow are not parallel to the nozzle axis. This results in a radial 
component of momentum which reduces the net axial thrust. A 
measure of this effect can be obtained by comparing the TDK 
one-dimensional kinetics solution to the two-dimensional 
kinetics solution. 

^sp(DIV) = I sp(ODK) “ I sp(TDK) 

The divergence loss can generally be minimized by carefully 
shaping the nozzle as discussed, for example, by Rao (1960) 
Lengthening the nozzle is also an alternative, though this 
may increase the nozzle weight to a point that it is not 
practical from mission considerations. 

The last loss mechanism is due to two types of viscous 
effects. The first is due to viscous drag on the nozzle walls. 
The second comes from flow blockage due to boundary layer 
growth, especially in the throat region. A measure of these 
combined effects can be defined by the difference in the TDK 
inviscid two-dimensional kinetics and boundary layer 
solutions 

^sp( VIS) = *sp(TDK) ~ I sp(BLM) ^ 

The shorter the nozzle, the less surface area for viscous drag 
and heat transfer to act on. Therefore, the viscous loss is 
minimized by reducing the nozzle’s wetted surface. The 
blockage effect can -be significant for small throat solar 
thermal nozzles. Its impact can be reduced by increasing the 
throat diameter. A more detailed study of this blockage effect 
will be the subject of future work. 

The overall nozzle performance loss is the sum of these 
individual losses. But, minimizing one of the individual 
losses may produce a corresponding increase in the loss from 
another mechanism. Therefore, maximizing overall nozzle 
efficiency is a process of compromise between each of the loss 
mechanisms. This paper describes the impact of geometric 
parameters on these loss mechanisms with the ultimate goal of 
developing guidelines for optimizing the total nozzle 
efficiency. 


2.3 Nozzle Design Constraints 

A baseline nozzle configuration and operating conditions 
were provided by Rocketdyne for this study. A sketch of the 
baseline nozzle with the geometric parameters for the TDK 
model of the baseline are shown in Fig. 3. Plenum chamber 
conditions were 50 psia and 6100°R (0.345 MPa and 3611 K). 
This temperature is the upper limit for a hafnium-carbide 



FIGURE 3: BASELINE NOZZLE. ALL DIMENSIONS 
NORMALIZED BY THROAT RADIUS, R*. 


absorber. The chamber pressure is limited by the yield stress 
of this material at the elevated temperature (Shoji, 1993, 
McClanahan, et al. 1994). These conditions produce an 
equilibrium molecular hydrogen mass fraction of 0.878. The 
nozzle throat area is 8. 43x10' 5 ft^ (7.83x 10'^ m^). The 
baseline area ratio is 100. This area ratio and chamber 
conditions produce approximately 1.0 lbf nominal thrust. In 
this study, the nozzle throat diameter, contraction area ratio, 
expansion area ratio and chamber conditions were all held 
constant. The four parameters that describe the throat contour, 
R u , Rj, 0 U , and 0,j. were varied independently. This 
combination of fixed and variable parameters allowed the 
convergent region cone angle and length to vary (with 0 U ), 
and allowed the divergent section contour (skewed parabola) to 
vary as well. The divergent section length was maintained 
constant at 0.1392 ft ( 0.0424 m). The boundary layer was 
allowed to transition from laminar to turbulent, with a criteria 
for transition to turbulence of Re > 400, with Reynold's 
Number based on momentum thickness. For all cases run in 
this study the boundary layer remained laminar. 


3 HYDROGEN DISSOCIATION RATE STUDY 

In a solar thermal rocket, liquid hydrogen is vaporized and 
partially dissociated. The recovery of kinetic energy through 
the exothermic recombination of hydrogen atoms increases 
the I S p of the nozzle. As discussed, nonequilibrium (finite 
rate) hydrogen kinetics can significantly affect this process. 
Assuming a pure hydrogen propellant and neglecting radiation 
effects, the exothermic recombination of hydrogen is 
governed by two reactions 

H + H + Mj <=>H 2 +Mj (5) 

where Mj is either an atomic (H) or molecular (H 2 ) third- 
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TABLE 1: HYDROGEN KINETIC REACTION RATE 
COEFFICIENTS (EQ. 6) 


Reference 

M { 

Aj 

N{ 

Bi 

Cohen & 
Westberg (1983) 





Low 

h 2 

0.5xl0 17 

-0.6 

0 

Nominal 


1.02xl0 17 

-0.6 

0 

High 


2.0xl0 17 

-0.6 

0 

Low 

H 

l.OxlO 15 

0 

0 

Nominal 


3.19xl0 15 

0 

0 

High 


l.OxlO 16 

0 

0 

NASP (1990) 

h 2 

1.83x!0 18 

-1.0 

0 


H 

7.3xl0 17 

-1.0 

0 

TDK Standard 

h 2 

2.56xl0 18 

-1.0 

0 


H 

1.6xl0 19 

-1.0 

0 

Baulch 

h 2 

2.98xl0 16 

-0.45 

-1.5 


H 

9.43xl0 16 

-0.34 

1.2 

Wetzel & 
Solomon (1993) 





Low 

h 2 

2.7xl0 16 

-0.58 

0 

Nominal 


9.7xl0 16 

-0.6 

0 

High 


2.0xl0 16 

0 

0 

Low 

H 

7.4xl0 16 

-0.63 

0 

Nominal 


9.1xlO i6 

-0.37 

0 

High 


6.0xl0 16 

0 

0 

FDNS 

h 2 

5.0x10* 5 

0 

0 


H 

5.0xl0 15 

0 

0 


Rate (&,•) Units = cm^/Mole 2 -sec 


body. The rate of each reaction can be written in the Arrhenius 
form 


( lOOOBj 'l 

kj=A|T N| e^ RT K (6) 

where T is temperature and R is the universal gas constant. 
Variation in experimental results have produced uncertainties 
in these reaction rates which are as much as two orders of 
magnitude (Wetzel and Solomon, 1993; Gerrish and Doughty, 
1993). This uncertainty may have a significant impact on the 
prediction of solar thermal rocket performance. The research 
literature offers a range of hydrogen third-body recombination 
rates in the temperature and pressure range of this study. These 
rates are summarized in Table 1 in the form of the constants for 
the Arrhenius rate of Eq. (6). 

Cohen and Westberg (1983) surveyed the literature and 
defined representative low, nominal and high rates. Research 
into NASP propulsion technologies resulted in a survey of 



Nominal Rate 


FIGURE 4: SENSITIVITY OF SPECIFIC IMPULSE TO 
REACTION RATE. 


hydrogen reaction rates and establishment of a standard set 
(Oldenborg, 1990). The TDK code (Nickerson, et al. 1993) 
and the NASA developed Finite Difference Navier-Stokes 
(FDNS) code (Chen, 1993) provide internal default rates. 
Rocketdyne provided baseline rates based on the work of 
Baulch, et al. (1972). Finally, Wetzel and Solomon (1993) 
studied the impact of rate uncertainty on the performance of 
high thrust (~ 10,000 lbf or ~ 45,000 N) nuclear thermal 
rockets. They also defined a set of high, nominal and low 
rates. The TDK code was run with each of these rates and the 
baseline nozzle under inviscid conditions. The resultant I cn 
values were calculated and are plotted in Fig. 4 as nominal 
rates (kxl0°). The rate uncertainty produces an uncertainty in 
I S p of approximately 50 seconds. Neglecting the Wetzel and 
Solomon high rates, the Baulch et al. (1972) rates produce an 
I s p which is approximately the average of the values 
calculated using the other rates. Therefore, these rates were 
chosen as the nominal rates for the rest of this investigation. 

The TDK code was next used to calculate the impact on the 
baseline nozzle I S p of increasing and decreasing the nominal 
rates by successive factors of 10. The baseline nozzle was 
also run with frozen and equilibrium chemistry specified. The 
analysis was made for both inviscid and viscous cases. The 
resultant curves are also shown in Fig. 4. These results 
illustrate that the nominal rates produce near frozen chemistry 
in the flow, indicating that most of the thermal energy used to 
dissociate the H 2 is not recovered in the nozzle, thereby 
reducing performance. Finite rate kinetics cause an I S p loss of 
approximately 170 sec (~ 13.8%) from the ideal equilibrium 
chemistry value. It should be noted that the magnitude of the 
kinetic losses is strongly dependent on the amount of 
dissociation in the chamber, and is therefore strongly 
dependent on chamber temperature. Lower chamber 
temperatures reduce kinetic effects significantly. The 
magnitude of performance lost due to viscous effects (210 sec, 
-18%) is comparable to the kinetics loss. 
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FIGURE 5: REPRESENTATIVE NOZZLE CONTOURS. 


4 EFFECTS OF GEOMETRIC PARAMETERS ON 
PERFORMANCE 

Parametric variation of the nozzle geometry was done with 
three goals in mind. First, the primary geometric parameters 
(i.e. length, throat diameter, area ratio) and inlet (plenum) 
conditions were held constant so that thrust levels for the 
various configurations were comparable. Second, 
recombination was to be maximized to increase I S p. Third, 
viscous effects should be minimized to increase Isp. 

Results obtained from studies of nuclear thermal rockets 
(Wetzel and Solomon, 1993; Gerrish and Doughty 1993) 
indicate that lengthening of the throat region encourages 
recombination. Therefore, the four parameters that define the 
throat contour, R u , R^, 0 U , and 0^, were varied independently. 
Figure 5 shows representative nozzle contours produced by the 
variation of these four throat parameters. The nozzle length 
was essentially maintained constant, except for cases of small 
0 U and large R u where the convergent section length is 
significantly longer than for the other nozzles. These longer 
nozzles do not maintain the original geometry, but are 
presented for completeness. 

A matrix of nozzle configurations were analyzed by the TDK 
code to determine the impact of geometry on I S p. The results 



Upstream Throat Angle, 0 U 


FIGURE 6a: PERFORMANCE TRENDS FOR VARIATION IN 

0 U - 



Downstream Throat Angle, e d 

FIGURE 6b: PERFORMANCE TRENDS FOR VARIATION IN 

0d- 


are shown in Figs. 6a - 6d. The ideal I S p of 1241 seconds was 
determined by a one-dimensional equilibrium (ODE) TDK 
analysis. The effect of each geometric parameter on the 
kinetic, divergence and viscous losses is expressed in terms of 
the AI S p definitions described in section 2.2. Improved 
performance is indicated by smaller values of AI S p. For the 
viscous calculations, both an adiabatic and constant 
temperature wall were analyzed. The constant temperature wall 
condition was set to 2460°R. This value was determined based 
on inclusion of radiation cooling through the exit plane and 
no liquid wall cooling. 

The I S p loss due to finite rate kinetics (ODE-ODK) is 
between 160 and 170 seconds. In these small nozzles this is 
primarily due to the relatively short stay times in flow regions 
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Downstream Throat Radius, R d 

FIGURE 6c: PERFORMANCE TRENDS FOR VARIATION IN 

Rd- 


where the temperature is decreasing and the pressure is still 
high enough to promote recombination. The loss is therefore 
minimized by increasing the throat length through decreasing 
the converging and diverging angles or by increasing the 
upstream and downstream radii. But, as shown in Figures 6a- 
6d, the kinetic loss is nearly insensitive to throat geometry. 
This result is not surprising since the reaction rates produce 
nearly frozen chemistry as shown in Fig. 4. 

This point is further illustrated in Fig. 7 where plots of 
atomic hydrogen mass fraction versus axial area ratio 
variation are shown for several representative geometries. 
These geometries were chosen to show the entire range of 
mass fraction variation found in the study. This figure 
indicates that the flow chemistry effectively freezes at 
approximately the same area ratio of 3:1 (and consequently 
close to the same pressure) for all nozzle geometries in this 
study. For the best case (R u = 10) only 1% of the mass is 
converted from atomic - to molecular hydrogen. 

Divergence losses due to two-dimensional flow (ODK-TDK) 
vary between 3 and 35 seconds. Both upstream geometry 
parameters have no impact on this loss, as expected. 

O 

Increasing the downstream diverging angle, 0,j, to 65 (see 
Fig. 6b) almost eliminates the divergence loss. Higher values 
of 0^ expand the flow quickly and allow more time for axial 
turning of the flow, while the lowest values of 0^ produce an 
almost conical nozzle, which is dominated by divergence 
losses. A lesser reduction in divergence loss is obtained by 
decreasing the downstream radius of curvature (see Fig. 6c). 
These effects somewhat compensate for the kinetic losses. 

Although modifications of the four geometric throat 
parameters have a relatively modest impact on the kinetic and 
divergence losses, the viscous losses can be significantly 
affected. In general, geometric changes which increase the 
wetted nozzle surface area increase the viscous losses. Viscous 
losses for the isothermal wall are typically larger than for the 
adiabatic wall. As shown in Fig. 6a, decreasing the upstream 
convergence angle produces a very long convergent section. 
This allows greater boundary layer growth and provides more 



Upstream Throat Radius, R 

FIGURE 6d: PERFORMANCE TRENDS FOR VARIATION IN 

Ru- 



FIGURE 7: REPRESENTATIVE AXIAL VARIATIONS IN 
MASS FRACTION OF H. 


o 

surface area for heat transfer losses. For 0 U > 45 , the 
convergence angle does not significantly affect the isothermal 
wall viscous loss. The adiabatic wall viscous loss decreases 
slightly with decreasing 0 U . Decreasing the downstream 
divergence angle, 0^, from 65° to 20° decreases the viscous 
losses by as much as 130 seconds of a total of 230 seconds 
(see Fig. 6b). This performance gain is offset slightly by the 
increase in divergence losses with decreasing 0^. 

Increasing the downstream throat radius of curvature, R^, 
decreases the viscous I S p loss by approximately 10 seconds. 
This gain is offset by the opposite trend in divergence loss. 
Increasing the upstream throat radius of curvature, R u , 
increases the viscous loss by approximately 10 seconds. But, 
this decrease in performance is offset by a comparable 
reduction in kinetic losses with increasing R u . 
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5 CONCLUSIONS 

Careful shaping of a nozzle within fixed constraints can 
potentially result in large gains in performance. This paper 
has dealt with the effect of nozzle geometry on three types of 
performance losses; finite rate dissociation-recombination 
kinetic losses, two-dimensional axisymmetric divergence 
losses, and compressible viscous boundary layer losses. 

The dissociation/recombination chemistry in these small 
nozzles was found to be nearly frozen. The lack of 
recombination caused an I S p loss of ~150 seconds when 
compared to equilibrium chemistry. Geometry variations had 
little effect on the chemistry. The effect may be greater and is 
worth investigating for larger nozzles. 

Divergence losses are well understood and documented in the 
literature. These losses accounted for no more than 30 seconds 
of I S p, and can be nearly eliminated through careful nozzle 
shaping. 

The boundary layer was the largest contributor to 
performance loss (AI S p ~ 200-300 seconds), but also showed 
the only significant reduction of losses with geometry 
variations. These losses were largely insensitive to upstream 
throat variations and downstream throat radius. The major 
gains were seen with variations in downstream throat angle, 
which largely determines the rate of expansion of the flow. 
The effect of the boundary layer was reduced when the 
downstream throat angle was low and the flow was allowed to 
expand more gradually. An I sp gain of -130 seconds was 
predicted by TDK when 0^ was reduced from the baseline value 
of 45° to 20°. 

There is some concern with the uncoupled boundary layer 
solver and its applicability to this problem. The assumption 
of an uncoupled inviscid core and viscous boundary layer 
solution is approximate, and generally deemed appropriate for 
'thin' boundary layers. The boundary layer in the small solar 
thermal nozzle occupies a large portion of the flow path. 
Further study needs to be done using more accurate methods or 
fully coupled solutions (i.e. CFD) to determine the magnitude 
of these trends. 

The numerical results of this study provide useful guidelines 
for the optimal design of low thrust solar thermal nozzles. All 
of these analyses should be verified with experimental data. A 
priority of future work should be placed on the design of 
appropriate experiments to validate these results. 
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